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Abstract: Many riparian (Tugai) forests growing along rivers in arid and hyper-arid regions of Central Asia 
are dominated by the Euphrates poplar (Populus euphratica). Besides generative reproduction, which is only 
possible upon flooding events and at a distance to the groundwater of less than 2 m, this phreatophytic 
tree species also reproduces vegetatively by forming clones that can cover land surface areas of several 
hectares. Along a gradient of groundwater distances, we investigated whether the fraction of clones in P. 
euphratica stands (1) increases with increasing distance to the water table; (2) is higher if supplied with 
water via river cut-offs; and (3) approaches 100% at a short distance to the groundwater, but at high salt 
concentrations in the upper soil layers, which would prevent germination and establishment of seedlings. 
AFLP (Amplified Fragment Length Polymorphism) analyses were conducted on leaf samples taken from 
mature P. euphratica trees growing at the fringes of the Taklimakan Desert in stands with different 
distances (2—12 m) to the groundwater at two plots at the middle and the lower reaches of the Tarim River 
and in a stand close to Ebinur Lake, Xinjiang, China. Genetic diversity was large among plots, but 
considerably smaller within plots. We found the highest genetic diversity (caused by regeneration from 
seeds) at plots that have a short distance to the groundwater or are supplied with additional water. There 
was no significant relationship between groundwater distance and clonal fraction. All investigated trees at 
the saline Ebinur Lake site belonged to one single clone. Our results demonstrate that the genetic pattern 
of this widespread species is not easily predictable even over small distances as it is a result of a complex 
interplay of stand history and dispersal of propagules (pollen, seeds, and vegetative diaspores) by wind 
and water. In conservation and restoration schemes, P. euphratica stands with a high genetic diversity and 
stands that grow at short distances to the water table and are regularly subjected to flooding (which favors 
generative over clonal reproduction) should be prioritized. 
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1 Introduction 


Every natural population, except for those that consist of one single clone, contains genetic 
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differences among their individuals. This also applies to most tree populations of forests. The 
genetic structures of those populations are produced by dispersal of pollen and seeds and thus are 
largely influenced by the mobility of these propagules. Besides generative reproduction, many 
tree species, however, also produce vegetative offspring. Vegetative reproduction of woody 
species, including alien and invasive taxa, can be found in forests and woodlands of Europe, 
South Africa, Asia, Australia, and North and South America, and plays an important role in the 
establishment and spread of populations (Ducci and Santi, 1997; Cremer, 2003; Douhovnikoff et 
al., 2005; Lorenzo et al., 2010; Milton and Dean, 2010; Vonlanthen et al., 2010; Thomas et al., 
2012). This reproductive strategy is particularly advantageous at sites with a harsh climate 
(Reisch et al., 2007; Eusemann et al., 2013), under unpredictable environmental conditions 
(Karrenberg et al., 2002) or after disturbance (Bellingham and Sparrow, 2000). At such sites, 
vegetative reproduction is often more efficient and successful than generative recruitment because 
clones have a higher chance to establish due to the constant supply with resources from the 
mother plant (Crawford, 2008). Among woody species, the creosote bush (Larrea tridentata) of 
the Mojave Desert in southwestern North America is known for its capability to form clones that 
reach a diameter of up to 22 m and an estimated age of almost 12x10? years (Vasek, 1980). In 
some cases, natural tree stands are even dominated by clonal structures. One of the most 
prominent examples for clonal tree colonies is the "Pando" clone of Populus tremuloides in 
western North America (Mock et al., 2008) with an area of 43 hm? and an estimated age of 
80x103 years. In general, vegetative reproduction is widely distributed in Populus (poplar) and 
Salix (willow) species of the willow family (Salicaceae) (Karrenberg et al., 2002; Reisch et al., 
2007; Stamati et al., 2007; Slavov et al., 2010; Eusemann et al., 2013; Tiedemann and Rood, 
2015). 

One of the habitat types favoring vegetative reproduction strategies is riparian forests along 
rivers in Central Asia. These forests are located in an extreme environment in close proximity to 
deserts and are often additionally subjected to high unpredictability of water supply, which 
depends on the amount of melting water from ice and snow in the adjoining high mountain ranges 
(Thomas et al., 2000). Hence, in these ecosystems, vegetative growth and the development of 
large stands of clonal trees provides distinct advantages for the establishment, distribution and 
interspecific competition of such species. Monospecific forests of Euphrates poplar (Populus 
euphratica and Populus diversifolia) along the Tarim River at the northern fringes of the 
Taklimakan Desert and close to the shore of Ebinur Lake in the arid northwestern part of China, 
forming extensive stands along the river banks and lake shores, constitute an ideal study system to 
analyze the genetic structuring and to address questions of clonal growth and gene flow in and 
among different stands of these trees. 

The studied species P. euphratica is a broad-leaved deciduous tree with a wide distribution in 
Central and Western Asia (Fang et al., 2015). In northwestern China, where it has its most 
extensive distribution, it is the predominant tree species of riparian (Tugai) forests along streams 
at the fringes of the Taklimakan Desert in southwestern Xinjiang (Wang et al., 1996). It also 
occurs near the shores of Ebinur Lake located in the southwestern part of the Junggar Basin in 
northeastern Xinjiang (Yu, 2008). According to the current United Nations Educational, Scientific 
and Cultural Organization (UNESCO) definition based on the ratio of average annual 
precipitation (P) to evapotranspiration (ET), the climate of these regions is hyper-arid (P/ET<0.03, 
fringes of the Taklimakan Desert; Xia et al., 1993; Thomas et al., 2000), or arid (0.03<P/ET<0.20, 
Ebinur Lake region; Yu, 2008). Under these climatic conditions, P. euphratica is an obligate 
phreatophyte, which—for survival and growth-relies on permanent access to the groundwater 
(Thomas, 2014). Thus, the P. euphratica-dominated Tugai forests in the arid and hyper-arid 
regions of Central Asia are unique riparian ecosystems that—due to a particular combination of 
climatic, geologic, hydrologic, and edaphic factors—hardly exist in other parts of the world. 
Deeper knowledge into the genetic structure of the forest stands and the rejuvenation capabilities 
of the trees is a decisive prerequisite for conserving these forests, which not only provide 
important ecosystem services to the local and regional population (Lang et al., 2015), but also 
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constitute a refuge for threatened animal species (Hai et al., 2006; Qiao et al., 2006). According to 
a previous study, these forests exhibit a large genetic diversity, but only small genetic 
differentiation among populations across distances of 9 degrees of latitude and 23 degrees of 
longitude in northwestern China (Wang et al., 2011). 

Populus euphratica forests establish when the water content of the surface soil, after an 
inundation by river or lake water, is high enough to allow germination of seeds and rapid 
formation of shoots and, most importantly, root systems. Generally, this happens after summer 
floods along river banks or lake shores with low salt concentrations (Thevs et al., 2008a). Upon 
continuous desiccation of the upper soil layers, the roots are forced to rapidly grow deeper into 
the soil to gain access to the groundwater or its capillary fringe to ensure survival and growth 
(Thevs et al., 2008b). Once being established, the trees can tolerate lowering of the groundwater 
level to a certain extent. Such a lowering of the water table might be caused by (1) changes in the 
course of the river (Thevs et al., 2008a); (2) withdrawal of large amounts of water from the river 
or the groundwater for irrigation of extensive agricultural fields, as is the case along the middle 
and lower reaches of the Tarim River at the northern and eastern fringe of the Taklimakan Desert 
in southwestern Xinjiang (Feng et al., 2005); or (3) sand accumulation upon wind-borne sand drift 
(Xia et al., 1993). During the latter process, tree crowns of P. euphratica reach distances of up to 
25 m from the water table (Gries et al., 2003). 

If a site is not flooded and the distance to the water table is larger than 2 m, no generative 
regeneration is possible because (1) seeds barely germinate on unflooded soils (Cao et al., 2012), 
and (2) the groundwater distance is too large to be bridged by the roots of the seedlings (Thevs et 
al., 2008a). Germination of P. euphratica is also hampered by increased salt concentrations of the 
substrate (Yiliminuer et al., 2015). Under such conditions, only vegetative regeneration is feasible 
by the formation of "suckers". These are shoots emerging from horizontally running underground 
shoot sections that are connected to the mother plant; this is a common feature of Populus species 
(Jeník, 1994). Such shoots are initiated in soil depths down to 0.6 m and reach distances of up to 
40.0 m from the parent tree (Wiehle et al., 2009). As a result of vegetative regeneration, large 
clones covering areas up to 121 hm? are formed (Vonlanthen et al., 2010). However, the trees' 
capability to regenerate vegetatively through suckers also decreases with increasing distance to 
the water table and is limited to a maximum groundwater distance of 8 m (Zerbe and Thevs, 
2011). Therefore, P. euphratica stands growing at a large distance to the groundwater often are 
relatively old and sparse, whereas dense poplar stands containing a large fraction of young trees 
are only found at sites with a relatively small distance to the water table or to the riverbed 
(Thomas et al., 2017). 

Some studies on the genetic structure and clonal growth of riparian P. euphratica forests and 
their relationship to the groundwater level have already been conducted in arid Northwest China. 
However, these studies were constricted to individual regions along the Tarim River (Eusemann et 
al., 2013; Petzold et al., 2013; Ling et al., 2015) and covered-if at all-a gradient of groundwater 
distances of up to 5 m only. Moreover, the age of the investigated trees has not been determined in 
these studies. Accordingly, Eusemann et al. (2013) pointed out that more extensive work is 
needed to provide evidence for relationships between clonality and the dynamics of the river and 
the groundwater. Thus, in our present study, we investigated the genetic structure of each three P. 
euphratica stands with known tree ages and different distances to the water table at the middle 
and the lower reaches of the Tarim River (with a distance of 360 km between these two plots) in 
southwestern Xinjiang and compared it to the genetic structure of a P. euphratica stand growing 
under saline conditions, but close to the groundwater at Ebinur Lake in the western part of 
northeastern Xinjiang (with a distance of 440-460 km to the plots at the Tarim River). The 
settings at the two Tarim River plots further allow assessing the impact of the actual and ancient 
river courses and the wind systems influencing gene flow among stands. 

In accordance with the life history of P. euphratica stands forming the riparian forests along 
rivers and lakes in arid regions of Central Asia, we formulate the following research hypotheses: 
(1) the fraction of clones in a given stand increases with increasing distance to the water table 
along a gradient from a distance to the groundwater at which generative regeneration is still 
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possible after flooding (2.0 m) to a groundwater distance that exceeds the limit of both generative 
and vegetative regeneration (12.0 m). However, in pondering the effects of the groundwater 
distance on the extent of clonality, we have to consider that the course of the Tarim River and the 
water flow through its cut-offs exhibited dynamic changes in the past (Thevs et al., 2008a; Yu et 
al., 2016). Therefore, we also hypothesize that (2) in stands with intermediate distances to the 
water table (5.0-8.0 m), the fraction of clones is larger in those stands that are additionally 
supplied with water via river cut-offs. This additional water supply should have promoted the 
spread of clonal structures upon a subsequent rise of the water table. Further, we expect that (3) 
clonal growth of P. euphratica predominates in a forest on a previous lake shore (Ebinur Lake) 
with a relatively short distance to the groundwater (1.5—3.0 m) because this forest has not been 
subjected to inundation for several decades and exhibits high salt concentrations in the upper soil 
layers. These high salt concentrations are due to a shrinking water body of the lake upon 
withdrawal of large amounts of water for irrigation of agricultural fields (Yu, 2008). 
Consequently, reproduction from seeds should not have been possible under those conditions for a 
longer period of time and the formation of new shoots should have been feasible only through 
vegetative recruitment. 


2 Materials and methods 


2.1 Study sites and sampling 


The study was conducted within the framework of the joint Sino-German project "Sustainable 
Management of River Oases along the Tarim River (SuMaRiO)" (2011-2016) at the middle 
reaches (Y, near to the village of Yingbazar) and the lower reaches (A, close to the village of 
Arghan) of the Tarim River, which runs along the northern and eastern fringe of the Tarim Basin 
in the southwestern part of the Xinjiang Uygur Autonomous Region, China (Table 1). 


Table 1 Location, distance to the groundwater, and variables of stand structure and tree morphology of the 
Populus euphratica stands at the study plots of Yingbazar (Y), Arghan (A) and Ebinur Lake 


YI Y2 Y3 Al A2 A3 Ebinur Lake 
Latitude 41°13.576'N 41°17.546'N  41°18.711'N 40°08.795'N  40°08.067'N_ 40°08.267'N__ 44°37.024'N 
Longitude 84°19.611'E 84°18.263'E 84°23.012'E 88°20.922'E 88°21.862'E  88°20.858'E  83°33.883'E 
Elevation (m) 930 931 929 829 829 830 280 
AMT (°C) 10.6! 10.6! 10.6! 11.72 11.72 11.72 6.84 
MAP (mm) 43! 43! 43! 153 153 153 1034 
Faa creer > 2500! > 2500! >2500! 2500-3000? 2500-3000? 2500-3000? 16624 
ERPE orsiance 2.0 7.5 12.0 5.2+EW 4.8 6.6 1.4-2.7 
Number or trees per 467 378 67 166 257 59 185 
hectare 
sso as 4 20 19 43 30 26 21 
Age óf ei trees 53,548.34 80.3+6.4°  —126.1420.3"  43.0+3.6° T1.545.8° 75.848.6° ~ 1124 
Height or sampled 16.5+0.9° 9.1+0.2° 6.9+0.4° 8.5+0.28 7.340.3° 6.0+0.2° 5.740.2 
trees (m) 
saris a 0.2740.02° 0.28+0.01°  0.62+0.03* —0.3140.02" —0.28+0.01"  0.33+0.02°  0.30+0.03 


Note: AMT, annual mean temperature; MAP, mean annual precipitation; DBH, diameter at breast height. EW, supply with "ecological 
water" in most years since 2010. Values that share the same lowercase letter do not differ significantly among plots at a given site (only 
one plot at Ebinur Lake). Mean+SE. !, Thevs et al. (2008a); 2, Gärtner et al. (2014); 3, Aishan et al. (2016); 4, Yu (2008); 5, in 2011 for 
Yingbazar, in 2012 for Arghan, and in 2015 for Ebinur. 

At each site, three P. euphratica stands (plots) were selected that differed in their distance to 
the groundwater or in water supply, which decreased from plot 1 to plot 3 at each site. In Arghan, 
two (Al and A2) out of the three stands exhibited approximately the same distance to the water 
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table (5 m), but in contrast to stand A2, stand Al is being supplied with additional "ecological 
water" in almost every year since 2000 for conservation and restoration of the poplar forests (Xu 
et al., 2006) through a nearby river cut-off, which, from time to time, may have also carried water 
in the years before. We assume that this water supply results in a temporary rise of the water table 
on plot Al. Under natural conditions, annual flooding at both plots generally occurs in late 
summer or autumn, that is, in the time of seed maturity (end of June through beginning of October; 
Wang et al., 1996). During the recent years and up to present, the study plots at both plots (Y and 
A) are uninfluenced by human activities or are influenced to a negligible extent as they are 
located in the protected Tarim Huyanglin Nature Reserve (Y1) or are too remote from settlements 
to be used for wood harvest or livestock grazing to a larger degree (Y2, Y3, and Al—A3). The 
location of all study plots and their distances to each other are provided in Figures S1—S3 of the 
Appendix. 

Each plot comprised a circular area with a radius of 50 m around a central P. euphratica tree 
(resulting in an area of 7854 m? per plot) and displayed a homogeneous stand structure. The 
distance to the groundwater was determined in spring (2012 in Yingbazar, and 2013 in Arghan) 
with a hand-driven soil auger (Edelman Sand; Eijkelkamp Soil & Water, Giesbeek, The 
Netherlands) at four (A1—A3) or seven spots (Y1); or with a drill hammer (Cobra mk1; Stitz 
GmbH, Gehrden, Germany), coupled to a hydraulic rod drawer (SZ—13A; Stitz GmbH), at one 
spot of plot Y2. In plot Y3, the distance to the water table was explored at one spot by a 
commercial company using the rotary drilling technique. 

The perimeter of each sampled tree was measured with a tape measure to calculate the diameter 
at breast height (DBH). The height of the trees was measured with an ultrasonic hypsometer 
(Vertex IV, Haglöf Co., Långsele, Sweden). Tree age was determined from stem increment cores 
or, in the case of some additionally selected trees (to increase the sample size), on the basis of the 
relationship between tree age and DBH according to Wang et al. (1996) and Thomas et al. (2017). 
No soil analyses were conducted for those two plots, but the influence on seed germination and 
vegetative reproduction by sodium and chloride in the upper soil is assumed to be low due to 
periodic flooding, which washes the salt out of the soil (plot Y1), or due to large distances of the 
soil surface to the capillary fringe (>4 m; plots Y2 and Y3 and all plots in Arghan), which 
prevents capillary flow of dissolved minerals from the groundwater to the upper soil layers. 

In addition to the study plots along the Tarim River, a fenced study plot in Ebinur Lake region 
with an area of 10,000 m? that had been used within the framework of a research project on forest 
development (Yu, 2008) was integrated in our study (Table 1 and Fig. S3 of the Appendix). On 
this plot, the water table was 1.4—2.7 m below the soil surface and the soil was characterized by 
high concentrations of sodium (12.1% at the soil surface) and chloride (390 mg/kg at a soil depth 
of 10 cm) as well as a high pH value (8.6 at 10 cm of soil depth) (Yu, 2008). Measurements of 
tree perimeter and height were carried out as described above. The average tree age of the stand 
was derived from Yu (2008). 

Leaves from approximately 20—40 trees per plot (Table 1) were sampled in early summer 2013 
and late summer 2014 in Arghan, late summer 2014 in Yingbazar, and June 2015 in Ebinur (at 
plot Y1, however, only 4 trees were sampled because the access to this plot was restricted for 
safety reasons upon administrational directive). From each tree, 20 leaves were sampled from the 
sun-exposed periphery of the crown. The leaves were air-dried and transported to the University 
of Trier (for analyses of nutrient concentrations in sub-samples) and afterwards to the 
Senckenberg German Entomological Institute for further processing and genetic analyses. 


2.2 DNA extraction 


Up to 30 mg of dried leaf material were pulverized with the TissueLyser LT (Qiagen, Hilden, 
Germany). Genomic DNA from the Arghan stands was extracted with the DNeasy Plant Mini Kit 
(Qiagen) and DNA from Ebinur and Yingbazar with the innuPREP Plant DNA Kit (Analytik Jena 
AG, Jena, Germany), using the OPT Lysis buffer. The extractions followed the manufacturer's 
instructions with the following minor modifications: to decrease possible inhibitor concentrations, 
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purified DNA was eluted first with 50 uL and, in a second step, with 100 uL elution buffer (in a 
separate tube). The 100 uL elution was used for further analyses. 


2.3 Amplified Fragment Length Polymorphism (AFLP) analysis 


The AFLP technique has been successfully employed to determine the clonal diversity and 
genetic structure of different plant species growing in harsh environments (de Witte et al., 2012). 
AFLP was made with fluorescence-labelled primers (6-FAM and HEX; Biomers, Ulm, Germany) 
and followed Vos et al. (1995). For each sample, genomic DNA concentration was determined 
with the QuantiFluor dsDNA System (Promega, Fitchburg, WI, USA) and standardized to 20 
ng/uL. For digestion, 100 ng of genomic DNA were used. The initial restriction-ligation step was 
performed for 16 h at 21°C using a thermocycler (Labnet, Edison, NJ, USA). Genomic DNA was 
digested simultaneously with restriction enzymes EcoRI and MselI. Double-stranded EcoRI and 
Mesel adapters were ligated to the sticky ends of the fragments. In the following two-step 
amplification, pre-selective primers with one selective base (EcoRI primer E, Msel primer M) and 
selective primers with two additional selective bases (EcoRI primer E+3, Msel primer M+3) were 
used. For the selective amplification, we used fluorescence-labelled primer combinations (Table 
S1 of the Appendix). 

AFLP analyses for the three Arghan stands were performed on a GE Healthcare 
96-capillary-sequencer (MegaBACE 1000) to separate fragments together with the internal size 
standard ET550-235 R (MegaBACE; GE Healthcare, Munich, Germany). The ramets from Ebinur 
Lake and Yingbazar were analyzed on a 3730XL DNA analyzer (Applied Biosystems, Waltham, 
MA, USA) together with the internal size standard 400HD (Applied Biosystems). Raw data were 
aligned with GENEMAPPER 5 (Applied Biosystems). All electropherograms were re-examined 
visually in order to check for possible misinterpretations of the automated GENEMAPPER 
analysis. Fragments were scored as present when the appropriate peak height exceeded the 
standard parameter-setting thresholds (blue: 60; green: 60). In the final binary matrix, the 
presence of a band was scored as '1', whereas the absence of a band was scored as '0'. The 
resulting matrix was used for all further statistical analyses. 

To check the reproducibility of the AFLP amplification and to estimate the error rate, we 
amplified and analyzed twice the eleven randomly selected samples. The error rate was calculated 
as the number of differences between the duplicated samples divided by the total number of 
fragments per profile (Bonin et al., 2004). 


2.4 Identification of possible clones 


In order to distinguish clones from closely related siblings and to calculate a threshold of 
similarity above which two AFLP profiles were considered to belong to one clone, we used the 
method of Douhovnikoff and Dodd (2003). We calculated the Jaccard similarity between the 11 
individuals analyzed twice and used for the estimation of the error rate; these were considered as 
"clones". As siblings, we randomly chose 25 individuals from Arghan and Yingbazar stands and 
determined the pairwise Jaccard indices in both cases. The threshold was calculated with the 
following formula: T=(SDcloneskmMeansibs+S DsipsxMeanclones)/(SDsibs+SDclones) (with T=threshold, 
SD=standard deviation, and sibs=siblings) (Schleuning et al., 2011). Hence, all pairs of 
individuals with Jaccard similarity equal or higher than T were assumed to be clones, all others 
not. 


2.5 Statistical analyses 


Within-stand genetic diversity was assessed using the percentage of polymorphic loci (%P pop), 
Shannon diversity index (H') based on all AFLP fragments of the dataset (Legendre and Legendre, 
1998) and Nei's gene diversity (Hn) (Nei, 1973). H' was calculated with the software POPGENE 
1.32 (Yeh and Boyle, 1997). Hy and the percentage of polymorphic loci were calculated for each 
stand with the R script AFLPDAT (Ehrich, 2006). The same script was used to calculate 
frequency down-weighted marker values (DW) (Schénswetter and Tribsch, 2005) with 
modifications given in Winkler et al. (2010). 
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The number of private and fixed private fragments (if present) per stand was extracted from the 
binary matrix. Private fragments are confined to a single stand. Fixed private fragments appear in 
every individual of a stand (Huck et al., 2009; Kramp et al., 2009). 

An analysis of molecular variance (AMOVA) (Excoffier et al., 1992) was conducted with the 
program ARLEQUIN 3.5.1.3 (Excoffier and Lischer, 2010). AMOVA values were calculated for 
the entire dataset and also for different genetic groups, in order to compare levels of 
among-stands differentiation between these groups. Significance levels of variance components 
were computed by a non-parametric permutation approach with 10,000 replicates. 

To examine similarities among and within stands, we subjected a matrix of pair-wise genetic 
distances among individuals of all populations to subsequent Principal Coordinates Analysis 
(PCoA). This analysis was performed using the Excel Add-in GENALEX 6.5b5 (Peakall and 
Smouse, 2012). 

From the binary matrix of the AFLP profiles, Nei's pair-wise genetic distances after Lynch and 
Milligan (1994) among all stands were calculated and relationships between them were 
constructed by applying the Neighbor Joining analysis (Saitou and Nei, 1987) using the software 
FAMD 1.31 (Schliiter and Harris, 2006). 

To study the underlying genetic structure in more detail, we used the multilocus assignment 
method of Pritchard et al. (2000) as implemented in the software STRUCTURE 2.3.4 (Hubisz et 
al., 2009). The program assumes that all individuals come from one or multiple unknown 
populations, each of them characterized by a set of allele frequencies at each locus. To perform 
inference, the method uses the Markov Chain Monte Carlo (MCMC) approach. STRUCTURE 
attempts to assign individuals to populations on the basis of their genotypes, while simultaneously 
estimating population allele frequencies (Pritchard et al., 2000). In our analysis, values of K (i.e., 
assumed genetic cluster) ranged from one to the number of the studied stands (i.e., seven). We 
used 10 replicates for every K with a length of burn-in period of 100,000 and a number of MCMC 
replicates after burn-in of 500,000. 


3 Results 


3.1 Stand structure and tree morphology 


According to our expectations, the stand density (number of trees per hectare) was higher and the 
average age of the sampled trees was lower on the plots with the shortest groundwater distance or 
conveyance of ecological water than on the plots with the largest groundwater distance at the 
study plots of Yingbazar and Arghan (Table 1; in Yingbazar site, statistically significant 
differences in tree age among the plots most probably were prevented by the small number of 
sampled trees on plot Y1). At both plots, the stems were shortest on the plots with the largest 
distance to the water table in Yingbazar and Arghan, and were also relatively short at the study 
site of Ebinur Lake. 


3.2 Population genetic parameters and structures 


The two selective AFLP primer combinations amplified a total of 269 fragments ranging from 63 
to 357 bp; 119 (i.e., 44.2%) of which were polymorphic fragments. Only well-defined and 
polymorphic fragments were selected for later clone detection (following the method of 
Douhovnikoff and Dodd (2003)). The technical error rate caused by scoring and labelling of the 
AFLP fragments was 2.66%. 

The analyzed stands of P. euphratica showed varying levels of genetic diversity (Table 2). The 
average number of polymorphic loci (#pol) was 107 in Arghan and 51 in Yingbazar site. The 
individuals from Ebinur Lake all showed an identical AFLP profile. The average gene diversity 
(Hn) was 0.115 (40.0066), and the average H' index was 0.134 (+0.079). The highest diversity 
was detected for Arghan site. The DW values showed very strong differences among the plots, 
ranging from 15.2 in Yingbazar to 137.1 in Arghan. 

PCoA based on genetic distances (Huff et al., 1993) assigned the three Arghan plots to three 
distinct genetic groups (coord 1: 24.9%, coord 2: 21.1%; Fig. 1a). However, the three Yingbazar 
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stands were genetically more mixed and showed no clear demarcations, but one of the individuals 
from plot Y2 was quite different from all the other trees analyzed from Yingbazar (coord 1: 
31.3%, coord 2: 8.2%; Fig. 1b). 


Table 2 Population codes, number of trees (N), number of polymorphic loci (# pol), Nei's gene diversity (Hn), 
Shannon index (H’), frequency-down-weighted marker values (DW) and number of (fixed) private fragments for 
the seven analyzed plots of P. euphratica 


Private fragment 


Population code N #pol Ay H' DW (fixed) 
Al 42 107.0 0.174 0.229 137.4 1 
A2 30 112.0 0.198 0.214 124.4 0 
A3 26 103.0 0.186 0.217 149.6 1 (1) 
Mean A 107.0 0.186 0.220 137.1 
SD 3.7 0.010 0.006 10.3 
Ebinur Lake 21 0.0 0.046 0.102 50.1 0 
Y1 4 95.0 0.090 0.061 4.1 1 
Y2 20 14.0 0.021 0.019 20.6 0 
Y3 19 44.0 0.094 0.096 21.1 0 
Mean Y 51.0 0.068 0.059 15.2 
SD 33.4 0.033 0.032 7.9 
Mean total 68.0 0.115 0.134 72.5 
SD total 44.0 0.066 0.079 57.8 
Note: Al, Al and A3 indicate Arghan stands. Y1, Y2 and Y3 indicate Yingbazar stands. SD, standard deviation. 
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Fig. 1 Associations among the individuals of the three Arghan (a) and three Yingbazar (b) stands of Populus 
euphratica revealed by Principal Coordinate Analysis (PCoA) performed on genetic distances calculated from 
AFLP data. Coord 1 and Coord 2 represent the first and second principal coordinates, respectively. 


A more detailed genetic structuring was revealed by the model-based Bayesian cluster analysis 
performed with STRUCTURE (Figs. 2a and b). The Evanno et al. (2005) method, implemented in 
the online program STRUCTURE HARVESTER (Earl and von Holdt, 2012), obtained the best 
resolution for K=7 (Figs. 2c and d), representing seven genetic groups found within all sampled 
individuals. Regarding the stands, Arghan 1 (A1) showed a rather mixed pattern including two 
genetic groups, one of them being partly reflected also in some individuals from Arghan 2 (A2) 
stand. Another noteworthy result was pointed out for Yingbazar 1 (Y1) stand, which entirely 
shared the same genetic cluster with Yingbazar 3 (Y3). 

An AMOVA test applied to the three stands from both Arghan and Yingbazar plots revealed 
similar levels of variation among and within stands. For Arghan, calculating the total variation 
gained 35.0% of variation among stands and 65.0% within stands. For the three Yingbazar stands, 
the variation was 35.7% among stands and 64.3% within stands. 
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3.3 Detection of clones within stands 


For the detection of clones, 119 well-defined fragments were selected (Douhovnikoff and Dodd, 

2003). Within the randomly chosen "siblings", pairwise similarity values were generated, ranging 

from 0.717 to 0.988 with a mean of 0.864 and a standard deviation of 0.086. Within the "clones", 
ie N (b) 
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Fig. 2 Genetic clustering analyses of P. euphratica genotypes. Results obtained from STRUCTURE analysis 
displaying the genetic groups found within all sampled stands. (a) geographic distribution of the P. euphratica 
stands analyzed in this study, including the individual affiliations to the seven genetic groups; (b) Neighbor 
Joining Tree based on Nei's pair-wise genetic distances; (c) result from STRUCTURE HARVESTER, displaying 
values of K (i.e., assumed genetic cluster)=7 being the best; (d) STRUCTURE bar plot of K=7. A1, A2, and A3 
indicate Arghan stands. Y1, Y2, and Y3 indicate Yingbazar stands. Ebi, Ebinur Lake stand. 


pairwise similarity values were generated, ranging from 0.978 to 1.000 with a mean of 0.991 and 
a standard deviation of 0.008. Hence, some pair-wise similarities among siblings were greater 
than pair-wise similarities among some of the clones. We established our boundary between 
clones and no-clones so that the tails of the standardized distributions were equal (see Materials 
and methods). This yielded a similarity threshold of 0.979. Based on this threshold, one (9.1%) 
clonal pairwise similarity value was lower than this threshold and two sibling pairwise similarity 
values (8.0%) were above this threshold (Fig. 3). 

Out of the 13,203 pairwise similarity values generated from the 163 individuals tested, we 
detected possible pairs of clones, using our similarity threshold criterion. Two pairs were found in 
A1, 5 in A3, 210 in Ebinur Lake (all pairwise similarities=1.000) and 154 in Y2; no signals for 
clones were obtained for A2 as well as Y1 and Y3 (Fig. S4 of the Appendix). 


4 Discussion 


4.1 Genetic structure and gene flow among stands 


AFLP analysis detected a generally low genetic diversity within stands, especially when 
compared to other studies on plants performed with this technique (Schönswetter and Tribsch, 
2005; Kramp et al., 2009). Thus, the genetic diversity and differentiation detected by our AFLP 
analysis differ strongly from the ones obtained by microsatellites (Wang et al., 2011; Eusemann et 
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al., 2013), in which the authors observed a high genetic diversity within stands. Furthermore, 
microsatellites revealed no genetic differentiation among stands, while AFLP detected gene pools 
of most of the stands analyzed in our study. Although, in principle, both methods are suitable for 
efficiently identifying polymorphisms (Miah et al., 2013), AFLP in this case might be the more 
appropriate marker system to reveal the genetic patterns across P. euphratica stands along 
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Fig. 3 Frequency histograms of pairwise Jaccard similarities for the 11 P. euphratica individuals that were 
analyzed twice and defined as "clones" (red bars) and 25 randomly chosen individuals treated as "siblings" (blue 
bars) to identify the similarity threshold (T) for distinguishing siblings and clones 


systems of river oases in northwestern China as it sensitively reflects the spatial distribution of the 
P. euphratica trees and stands. 

Hence, gene flow among different stands of P. euphratica growing at the same regional site 
might not be as high as previously assumed (Wang et al., 2011; Eusemann et al., 2013). In 
consequence, even stands separated by only some few kilometers can have their own genetic 
make-up, so that inter-stand gene exchange apparently is rare. Nevertheless, rare events of gene 
flow among local stands can occur. On the basis of STRUCTURE and PCoA, this is obvious from 
the existence of one 117-year-old tree growing in the A2 stand, which could represent an 
offspring between one parental individual from either Al or A3 and one parental individual from 
A2. This demonstrates successful pollen dispersal from either Al or A3 to A2 and germination of 
the resulting seed at a time when the distance to the groundwater and to a water-carrying part of 
the river was much shorter. Pollen can be easily transferred among the stands of a given site at the 
Tarim River as days with strong winds are frequent in the Tarim floodplain; those winds come not 
only from the northeast but also from the northwest and can even intensify into dust storms, 
especially in the eastern part of the Tarim Basin (Hai et al., 2006). Another possibility could be a 
possible hybrid between P. euphratica and P. prinosa, a species also distributed at the upper 
reaches of the Tarim River. 

Some few additional trees of the Al and A2 stands show varying amounts of the respective 
other gene pool and hence underline the successful reproduction of such inter-stand mixed 
individuals in their respective stands. One tree has signals of gene flow from A3 to Al in the 
PCoA analysis, but the signal is weak. Apparently, this inter-stand gene flow is not sufficiently 
large to level out the stand-specific genetic make-up in most cases. Exceptions might be Y1 and 
Y3, which showed no stand-specific genetic make-up, neither for the STRUCTURE nor for the 
PCoA analysis, in which trees from both stands appear genetically and completely mixed. The 
genetic similarity of the stands Y1 and Y3 could be due either to a constant gene flow between 
them or to the comparable small sampling size (four individuals) obtained from Y1. A possible 
reason for this mixture might be the transport of diaspores through river cut-offs that had 
interconnected the plots in former decades. Such a water-driven translocation of propagules may 
have also been the cause of the establishment of two Y2 ramets with an age of 60 and 80 years 
that, in most probability, represents one separate clone, but is genetically similar to the Y1 and Y3 
individuals. Remainders of cut-offs are still visible, but are now dry at least since 2011, possibly 
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due to changes in the local water management (including withdrawal of water for the irrigation of 
agricultural fields; Yimit et al., 2006). Altogether, our study demonstrates that the genetic 
structure of P. euphratica stands can exhibit a complex pattern as a combined result of 
pollen-mediated gene flow and clonal reproduction. A similar phenomenon has been observed in 
stands of P. trichocarpa in northwestern North America. There, sub-populations have been 
identified within a relatively small stand in spite of extensive pollen-mediated gene flow; this 
finding was explained by assuming different seedling cohorts produced by relatively few mothers 
(Slavov et al., 2010). 


4.2 Clonality of stands and their distance to the groundwater 


Our findings of strong superiority of generative over vegetative reproduction in the stands Al and 
Y1 with high groundwater table or additional water supply support our first hypothesis and the 
statements by Thevs et al. (2008b) that generative reproduction of P. euphratica depends on flood 
events and on a short distance to the groundwater. However, the low percentage of vegetatively 
reproducing individuals detected in A3 and Y3 is contrary to our expectation of a straightforward 
increase in the share of vegetative reproduction with an increase in the distance to the 
groundwater. This result parallels findings obtained from other Salicaceae species (Salix herbacea, 
Salix lanata, and Salix lapponum), which showed that even under the harshest conditions along an 
environmental gradient, sexual reproduction still prevails over clonal growth despite of the 
species’ high capacity for vegetative recruitment (Stamati et al., 2007). On the other hand, our 
data support the assumption by Wiehle et al. (2009) that a rather large distance to the water table 
reduces the fitness of clonally growing individuals, hence selecting in favor of generatively grown 
trees. Likewise, clones of creosotebush (Larrea tridentata) in the Mojave Desert of southwestern 
North America were larger and more common on substrate with beneficial soil water conditions, 
but rare and less dense at drier sites with a temporarily variable soil water regime (McAuliffe et 
al., 2007). At our study plots A3 and Y3 with a large distance to the groundwater, the small 
fraction of clones may be caused by a higher mortality of ramets due to increased hydraulic 
resistance along the paths of water flow from the groundwater to the tree crowns (Gries et al., 
2003; Thomas et al., 2008) or due to asymmetric competition for water among the ramets 
(Callaghan et al., 1992). Unexpectedly, most of the trees of Y2, a stand with an intermediate 
groundwater distance across the investigated gradient at that site, apparently belonged to a single 
clone well differentiated from Y1 and Y3. In spite of the relatively large distance to the 
groundwater (7.5 m), the Y2 stand exhibited a high tree density and a large fraction of young 
trees (46% of all trees on that plot belonged to the age class of 0-20 years; Thomas et al., 2017). 
This indicates that the water table had decreased within a relatively short period of time. 

With regard to the stands growing at intermediate groundwater distances, the picture obtained 
by our results is more differentiated than had been expected. According to general theory, our 
second hypothesis predicted a strong favoring of vegetative growth at intermediate distances to 
the groundwater because successful establishment of trees from germinating seeds is little likely 
under such conditions, but access to the groundwater is easy for well-established trees and water 
supply sufficiently secures that the risk of die-offs is rather limited for interconnected clones. 
However, our data follow this assumption only for Y2, but not for A2. The Tarim River is a 
highly dynamic system changing within decades considerably (Thevs et al., 2008a; Yu et al., 
2016) and the use of its water for irrigation purposes has increased tremendously over the past 
few decades (Feng et al., 2005; Xu et al., 2006). Thus, most likely, the water table on plot A2 was 
still high when most of this stand was established, but dropped considerably during the 
subsequent decades. This is even more plausible when taking into consideration that the stands in 
Arghan, located at the lower reaches of the Tarim River (as opposed to Yingbazar at the middle 
reaches), are much more affected by the increasing consumption of water for irrigation of arable 
fields (Peng et al., 2014). However, we also must not underestimate the influence of stochasticity 
in these extreme habitats so that occasional groundwater fluctuations and erratic inundations 
might be responsible for the structure of the stands in Arghan. It should be noted that the stands at 
plots with an intermediate groundwater depth have to be handled as a heterogeneous group with 
potentially large differences in their genetic stand structure. 
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According to our third hypothesis, our assumption of a clear favoring of clonal growth on salty 
soils with a groundwater depth (in parts, <2 m) that otherwise is well suitable for generative 
reproduction (Wang et al., 1996) is strongly supported by our finding of the highest extent of 
clonality in Ebinur Lake. 


5 Conclusions 


Our findings may have important consequences for the future conservation planning, and genetic 
analyses should be included for an optimization. As safeguarding genetic diversity is crucially 
important, the ecologically best secured stands at places with high groundwater level and frequent 
inundations should be given the highest conservation priority. These stands seem to be mostly 
composed of trees resulting from generative reproduction and, thus, are the result of genetic 
intermixture and hence uniqueness. Flooding events have already been proven successful in 
initiating the restoration of riparian forests at the lower reaches of the Tarim River. In the case of 
intermediate groundwater levels, genetic analyses are pivotal for the selection of the most relevant 
stands as only the ones representing genetically diverse stands with a high percentage of trees 
resulting from sexual reproduction may withstand further reduction in water supply, even if the 
shortage is just temporary. In case of stands at places with a large distance to the groundwater and 
unlikely inundation events, a clear focus has to be given to genetically diverse stands as stands 
mostly composed of clones might have a very poor survival perspective. For these reasons, we 
also recommend adapting the regional water management to safeguard the genetically most 
relevant stands. Against the background of limited water resources in the entire region, however, 
the supply with water for conservation purposes should concentrate on plots that have a short 
distance to the groundwater and can be flooded at least once a year to ensure the long-term 
survival of genetically diverse riparian poplar forests. 
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Fig. SI Yingbazar study site with the three study plots Y1, Y2, and Y3. The current Tarim River course can be 
seen at the lower part of the picture; at the upper part of the picture, old cut-offs between the study plots Y2 and 
Y3 are visible. These were once filled with water; but after 1950/1960, the water of the Tarim River was 
over-exploited due to the large increase in number and size of cotton farms (these farms are scattered throughout 
the entire study site) resulting in a drying-up of these cut-offs (Thevs et al., 2008b); thus, the water course 
between the study plots Y2 and Y3 disappeared. 


es 


Fig. S2 Arghan study site with the three study plots Al, A2 and A3. The picture shows the desert-like character 
of this site. The Tarim-River (displayed on the right side of the picture) and the cut-off above plot A1 are filled 
with water during ecological water diversion only (A1 is not inundated during ecological water diversion, but the 
groundwater level rises). Plots A2 and A3 do not benefit from ecological water diversion as they are too far away 
from the river and the cut-offs. 
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Fig. S3 Ebinur study site close to Ebinur Lake (visible on the upper left side of the picture) 
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Table S1 Adapter, primer, and fluorescent primer combinations used for the AFLP analysis of Populus 
euphratica 


Adapter Sequence 
EcoRI 5'-CTCGTAGACTGCGTACC-3'; 3'-CATCTGACGCATGGTTAA-S' 
Msel 5'-GACGATGAGTCCTGAG-3'; 3'-TACTCAGGACTCAT-S' 
E+A GACTGCGTACCAATTCA 
M+C GATGAGTCCTGAGTAAC 
sa SS 
E34 GACTGCGTACCAATTCAAT 
E45 GACTGCGTACCAATTCATG 
M54 GATGAGTCCTGAGTAACCT 
M61 GATGAGTCCTGAGTAACTG 
m a yS 
6-FAM E34/M54 (blue) 
HEX E45/M61 (green) 


Fig. S4 List of possible clone pairs. The Jaccard similarity indices >0.979 are marked in red. Possible clone 
pairs are listed for the Arghan (a), Ebinur (b), and Yingbazar (c) stands. The three-digit number is the 
identification of the individual Populus euphratica ramets sampled. 


